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By the Bethe ansatz method we study the energy dispersion of a particle interacting by a local 
interaction with fermions (or hard core bosons) of equal mass in a one dimensional lattice. We 
focus on the period of the Bloch oscillations which turns out to be related to the Fermi wavevector 
of the Fermi sea and in particular on how this dispersion emerges as a collective effect in the 
thermodynamic limit. We show by symmetry that the dispersion is temperature independent for a 
half- filled system. We also discuss the adiabatic coherent collective response of the particle to an 
applied field. 
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Prototype integrable quantum many-body problems 
provide valuable insight into the physics of correlated 
electronic and magnetic systems. With the development 
of novel experimental systems, as quasi-one dimensional 
quantum magnets [T] and cold atom systems [2], it be- 
came possible to tailor make and experimentally study 
these systems. At the same time because of their non- 
generic character they exhibit unconventional dynamic 
and transport properties of academic and potentially 
technological interest. 

One of the simplest models that was studied early on 
[3] , is that of a particle interacting with spinless fermions 
in a one dimensional lattice, where the mass of the ex- 
tra particle is the same as the mass of the fermions. In 
this model and its lattice version, it was shown that, 
despite the particle-fermion interaction, the transport 
of the "heavy" particle remains ballistic at all temper- 
atures, namely its mobility diverges [H [5]. The spin- 
less fermion or hard core boson bath corresponds to the 
Tonks- Girardeau gas, a one dimensional boson model 
with very strong local repulsion. 

More recently [IHl], motivated by cold atom physics 
experiments, the dispersion and period of Bloch oscilla- 
tions of the heavy particle when acted by a constant field 
in a continuous system were also studied and shown to 
be related to the Fermi wavevector kp of the fermionic 
sea. The first issue we study in this work is the compe- 
tition between the period of oscillations imposed by the 
lattice in the tight binding version of the model and that 
imposed by the Fermi wavevector. 

The lattice model corresponds to a particular sector of 
the ID Hubbard model [5], a prototype integrable model 
by the Bethe ansatz method [H]. The analytical solution 
allows us to characterize the low energy spectrum and 
also study the adiabatic response of the particle to an 
applied field. 



The system is described by the Hamiltonian, 

H = -thY,{^"^A+idi + h.c.)^tY,{c\+^ci + h.c.) 
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where c/(c|) are annihilation (creation) operators for N 
spinless fermions and di {dj ) of the extra particle on an L 
site chain with periodic boundary conditions. The inter- 
action comes only through the on-site repulsion U > 0. 
If the particle is charged and depends linearly on time, 
(j) — —£t, then a constant electric field £ — ~d<j)/dt acts 
on it. For a vanishing field (it could also be gravitational 
PI) the system should follow the adiabatic ground state 
that we will now map. 
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FIG. 1: Low energy spectrum for L = 10, u — 0.5, th = t, 
along with the reconsruction of selected branches by the Bethe 
ansatz solution. 

For illustration of the typical finite size lattice spectra. 
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we show in Fig. [T] the evolution of the low lying energy 
levels for L = 10, half- filling and in the k = Q symmetry 
momentum subspace as a function oi (f). It is clear that 
by scanning we recover at (/> = 21:111 L (n = 0, L — 1) 
the spectra in the succesive fc— subspaces with ^ = 0. 
We will now show how to reconstruct selected branches 
of the dispersion and study their finite size effects using 
the Bethe ansatz solution. 

The Bethe ansatz wavefunctions, in the presence of flux 
0, are characterized hy M — N quantum numbers kj 
for TV fermions plus the one extra particle, given by the 
following equations (we take = t), 
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+ -0(sinfc, -A), j = l,...,Af (1) 



6{p) = -2ia.\r^(p/u), u = U/At 
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The total energy and momentum are given by the sum 
of quasi-energies and momenta. 
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^6(fc,) = -2t^cos%, K = Y,k,. 



(4) 



Every state is characterized by a set of half-odd integers 
Ij and the (half-odd) integer J for (even) odd number 
of fermions. We will consider the case of odd number of 
fermions, M even, that is also equivalent to that of hard 
core bosons in a one dimensional system. 

We can reconstruct the low energy states by consider- 
ing the following branches; 

(a) place the I'^s at the successive values —M/2 + 1/2 < 
Ij < +M/2 — 1/2 and scan A between ±cxi. At this 
configuration, for u — >■ and A = the phase shift term 
in eq.Q takes the value +tt / L{—tt / L) for kj < 0{kj > 0). 
Thus for one particle kj = while the rest of fermions 
fill uniformly a Fermi sea, between —{N — l)/2 < kj < 
+{N — l)/2. The k-occupation corresponds to that of 
independent particles. 

For a finite u, by varying A every kj shifts by 
-7r/L(+7r/i) for A -cx3(-|-oo) so that -2Mit/L < 
K < +2MTr/L. In the thermodynamic limit —kp < 
K < +kp with kp — ttN/L. It is clear that the spec- 
trum is scanned by finding the A solutions of eqs.Q as 
a function of </> for J = or equivalently for = as a 
function of the quantum number J. 

(b) shift each Ij of branch (a) hy +1. Now by varying 
-oo < A < -l-oo, ttM/L < K < SttM/L which tends in 
the thermodynamic limit to kp < K < Skp- 



(c) shift any Ij of branch (a) by L and scan A. 

(d) shift Im of branch (a) by -1-1 and scan A. 

The adiabatically evolving ground state is obtained by 
the successive shift of Ij Ij + 1 up to L times when 
the k'jS are shifted by the trivial phase 2tt. Branch (d) 
corresponds to a state with an " electron-hole" excitation 
that is not adiabatically evolving from the ground state. 

Next, in Fig. [2] we show the finite size dependence of 
the adiabatially evolving ground state for a series of lat- 
tices at half-filling {kp = 7r/2). It is now clear that in the 
thermodynamic limit the period becomes 2kp. In order 




FIG. 2: Finite size scaling of adiabatically evolving ground 
state for th = t, u = 0.5 and half-filling. 



kj ~ fcO + ^6i(sin fc° - A), J = 1, . . . , M. 



(5) 



to analyze the finite size effects and obtain the particle 
dispersion we take to 0{1/L), 

1 

-( 
L 

We can now see how the dispersion becomes symmetric 
around kp in the thermodynamic limit by evaluating the 
difference in energy between the state A' = and K = 
2kp. By shifting fci = {2tt/L){-M/2 + 1/2) to ky = 
(27r/L)(M/2 + 1/2) and taking A = 0, we obtain, 
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K2k^ ~Ko^ —M + 0{l/L) ~ 2kp (6) 
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E2kF - Eo = +2icos(fci) - 2tcos{ki>) 
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sin fci?. (7) 



So the finite size effects decrease as 1/L. Next we can 
obtain an analytical solution of the dispersion of branch 
(a) and by symmetry (b) in the thermodynamic limit 
including corrections 0{1/L) from eqs. ( [3ppl ), 

M 

K = '^kj Kq + SK 
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;sinfcj"-A) (9) 



Replacing sums by integrals for the branch (a) we find 
for the momentum and energy dependendent terms on A 
and thus 6, 
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sin fc6'(sin fc — A). 



(10) 



Scanning — oo < A < +oo we obtain the dispersion shown 
in Fig. [2]for L ^ +oo. 

If we interpret SE, SK as the dispersion of the correla- 
tion energy and momentum we can deduce their temper- 
ature dependence by inserting in eqs. (8[9) a Fermi-Dirac 
thermal distribution for the momenta fc^. 
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Now, by electron-hole symmetry we see that at half- 
filling, fi = 0, the thermal distribution factors cancel so 
that the dispersion becomes temperature independent. 
This is due to the particular form of scattering phase 
shifts in this integrable model. 

The next issue we want to discuss is the response of the 
particle to a constant external field £ created by a time 
dependent fiux </> = —St. As is evident from Fig. [T]the 
collective ground state evolves through a maze of excited 
"electron- hole" states by level crossing. The coherent 
evolution throught the energy spectrum is protected by 
the macroscopic number of conservation laws character- 
istic of an integrable system that suppress level repulsion 
and also result to Poisson statistics. There is a coher- 
ent drag of the Fermi sea by the "heavy" particle. We 
should note that a similar level crossing trajectory also 
characterizes the time evolution of every excited state. 
We can plausibly argue that this is a generic behavior of 
integrable quantum many-body systems. 

To put this evolution into evidence we show in Fig. [3] 
the long time adiabatic energy evolution of an L = 10 
system as a function of phase </) = —£t in units of tt for 
£ — >■ 0. When the momentum of each particle is displaced 
by 27r so that (fi — 2ttM we have an identical state (in 
the example at (p/ir = 2M, M = 6). Thus the periodicity 
of the coherent collective motion becomes macroscopic. 
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FIG. 3; Time evolution for L 
th = 0.5t. 
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0.5 at th = t, and 



In the adiabatic limit there is reversible pumping of a 
macroscopic amount of energy into the system that it 
could experimentally be observed in a mesoscopic system. 
In contrast, in the nonintegrable case e.g. th — 0.5t, the 
periodicity is simply 2t: for a finite size system and it 
becomes 2k p in the thermodynamic limit. This behavior 
is generic for any interaction U and hopping difference 
tfi 7^ t. Note also that in the continuous system, hard core 
bosons in ID, the driving leads to an indefinite increase 
of energy as there is no = 2'kM periodicity. 

As shown in Fig. [3j for a finite field e.g. £ = 0.1, in the 
integrable case the energy deviates from the adiabatic - 
nonlinear in energy - trajectory, in sharp contrast to the 
nonadiabatic case where it spreads diffusively upwards 
from the ground state. It is clear that in the integrable 
system the departure from the adiabatic evolution is by 
mixing with highly excited states while in the noninte- 
grable case by tunneling between level repelled nearest 
neighbor states. For long times, the integrable system 
rapidly evolves through the whole energy spectrum to 
thermalize at the infinite temperature energy limit (mean 
value of energy spectrum) in a nonmonotonic way while 
in the nonintegrable system diffuses monotonically to the 
same mean energy (not shown) . Of course this unconven- 
tional thermalization is irrelevant in the thermodynamic 
limit but it might be observable in a finite size system. 

To quantify the spreading of the wavefunction during 
the evolution, we show in the Fig|4]A(0) defined by (TU], 

A(^) *(</,) I (H((/))- < H(0) >f\^{<t>) > (12) 

Here we also observe that the spreading is highly irregu- 
lar in the integrable case while it is almost linear in time, 
diffusive, in the nonintegrable one. 

The main issues emerging from these results are, first, 
whether it is possible to observe experimentally the adi- 
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ing to what extent this coherent behavior is affected by 
nonadiabatic effects, for at least some modes of driving. 
It is experimentally interesting to tune the boson-boson 
interaction from the weak (nonintegrable) to the strong 
Tonks- Girardeau (integrable) limit and observe the re- 
sulting dynamics. Or alternatively vary the "heavy" par- 
ticle mass. 

We conclude that this prototype integrable model ex- 
hibits an unconventional collective nonlinear evolution 
when driven by a constant field that motivates the ex- 
ploration of other similar integrable models. 

It is a pleasure to thank D. Gangardt and A. Kolovsky 
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032980-2 NOVMAG project. 



FIG. 4: Time evolution of wavefunction spread A((;/>) for L = 
10, u = 0.5 at th ^t,th = 0.5t and £ = 0.1. 

abatic dispersion of the driven particle in a macroscopic 
system, that is in fact qualitatively similar in the inte- 
grable and nonintegrable case. In this context we should 
mention that recent works argue on the absence of adi- 
abatic limit in low dimensional gapless systems [11 in 
macroscopic systems. Still, despite the spreading of the 
wavefunction, it might still be possible to observe the 
characteristic energy-momentum dependence for at least 
a few Bloch oscillations and the remnants of the singular 
integrable behavior. 

Second, in a finite size - mesoscopic - system, there is a 
striking difference between the nonlinear evolution of an 
integrable and a nonintegrable system. It is worth study- 
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